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ABSTRACT 



We examine the prospect to measure the Higgs boson mass using the recently introduced 
kinematic variable, the M^-Assisted On-Shell (MAOS) momentum, that provides a 
systematic approximation to the invisible neutrino momenta in dileptonic decays of 
jy-boson pair. For this purpose, we introduce a modified version of the MAOS 
momentum, that is applicable even when one or both of the jy-bosons from the Higgs 
decay are in off-shell. It is demonstrated that the MAOS Higgs mass distribution, 
constructed with the MAOS neutrino momenta, shows a clear peak at the true Higgs 
boson mass when an event cut selecting higher value of Mt2 is employed. We perform the 
likelihood analysis for this MAOS mass distribution to determine the Higgs boson mass, 
and find it can improve the accuracy of the Higgs mass measurement. Our results 
indicate that the MAOS Higgs mass can be useful also for the discovery or exclusion of 
the Higgs boson in certain mass range. 
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1 Introduction 



Hunting down the Higgs boson, the last ingredient of the Standard Model (SM), is one 
of the most important tasks of the LHC [1,2]. The LEP experiments established a lower 
bound of 114.4 GeV, at 95% confidence level, on the mass of the SM Higgs boson [3]. 
On the other hand, the electroweak precision data points towards a relatively light SM 
Higgs boson with run ~ 185 GeV at 95% confidence level [4]. Recently, Tevatron data 
have excluded the SM Higgs mass in the range 160 GeV < m# < 170 GeV again at 95% 
confidence level [5]. 

The strategy for the Higgs boson search depends on its decay pattern. The SM 
Higgs boson lighter than about 180 GeV mainly decays into the b quarks and W bosons. 
At hadron collider such as the Tevatron or the LHC, Higgs boson search through the b 
quark channel appears to be very difficult due to overwhelming QCD backgrounds. In this 
respect, the Higgs decay H — > WW — > Ivl'v' with I , /' = e , \i may provide the best search 
channel for the SM Higgs boson in the mass range 135 GeV < m# < 180 GeV. Even for 
a heavier Higgs boson with m# > 2Mz, this channel gets benefit from a larger branching 
ratio compared to the decay H — > ZZ — > 4 1. A drawback is that this channel involves two 
invisible neutrinos, making it impossible to reconstruct the Higgs boson mass directly. One 
then has to rely on a Higgs-induced excess in the distribution of certain observables, and 
this procedure typically requires an accurate estimate of the background contributions [6]. 

Recently, a new collider variable, the M^-Assisted On-Shell (MAOS) momentum, has 
been introduced [7] to approximate the invisible particle momenta in the process XX' — > 
VxV'x'i where X and X' denote pair-produced mother particles, V and V represent the 
visible particles (one or more particles for each of them) produced by the decays of X and 
X', respectively, and x an d x' are t ne invisible particles having the same mass. In this paper, 
we examine the possibility to determine the Higgs boson mass using the MAOS momenta 
of neutrinos in the process H — > WW — > Ivl'v' . An interesting feature of this approach is 
that one can use the collider variable Mj- 2 [8] for event selection, which enhances both the 
signal to background ratio and the efficiency of the MAOS momentum approximation to the 
true neutrino momentuirfl As we will see, the MAOS Higgs mass distribution, constructed 
with the MAOS neutrino momenta under a suitable M?2 cut, shows a clear peak over the 
background at the true Higgs boson mass. One can then determine the true Higgs boson 
mass by performing the likelihood fit to the MAOS Higgs mass distribution. The results of 
our likelihood analysis indicate that the precision can be significantly improved when the 
MAOS Higgs momentum is used for the Higgs mass determination, possibly combined with 
the kinematic variables considered in [1,10,11]. One can do a similar likelihood analysis for 
the Higgs boson discovery or exclusion, which might enhance the significance of the result 
for certain range of the Higgs boson mass. 

*It has been pointed out recently that Mti can be used also for the selection of new physics events [9]. 
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2 M^2-assisted on-shell (MAOS) momentum 



In this section, we discuss some features of the MAOS momentum for the dileptonic Higgs 
boson decay: 

W(p + k)W(q + l) -> l{p)u{k) l'(q) u'(l). (1) 

For the sake of discussion, we decompose the final state momenta into the transverse and 
longitudinal parts as follows: 



= (VIptI 2 +pLpt,pl), W=(J\k T \> + k? L7 k T ,k L 



t = WM 2 + ql,qT,q L ), l" = ( V I ^ 1 2 + l l , l T,h), (2) 

where we have neglected the masses of the charged leptons and neutrinos. For a Higgs boson 
mass tuh > 2M^, the two W bosons are in on-shell. On the other hand, if m# < 2M^, 
one or both of W bosons should be in off-shell. Regardless of whether the W bosons are in 
on-shell or not, one can construct the event-by-event variable Mt2 [8] which is given bjfl 



Mt2 = rnin 



max {M^\M^}] , (3) 



where pV denotes the missing transverse momentum carried by neutrinos, and are 
the transverse masses of the two decaying W bosons, which are given by 



^ 2 = 2(|p T ||k T | -pr-kr), (A4 2) ) 2 = 2(|q T ||l T | -q T -l T ). (4) 



If the Higgs boson is produced without having a sizable transverse momentum, so 
that the W-pair transverse momentum p^ w ~ 0, the missing transverse momentum is 
(approximately) given by 

i>T = - (pt + qr) • (5) 
In this case, M?2 is simply given by [12] 

M| 2 = (M^) 2 = (MP ) 2 = 2(|p r ||q T |+p T -q T ). (6) 

In fact, for both the signal W bosons from H —* WW and the background W bosons 
from qq — > WW, initial state radiations (ISR) typically give§ \pt W \ = 10 ~ 40 GeV, 
which would be non-negligible compared to the typical lepton transverse momenta in the 
final state. Still the above expression of Mti obtained in the limit p^ w = shows some 
qualitative feature of Mt2 for non-negligible Pj^ w , e.g. the correlation between Mt2 and 
the transverse opening angle A<3>^ between p^ and q^- In the following, ISR effects will be 
fully included in the numerical evaluation of M?2 and the MAOS momentum. 



TFor recent applications of Mt2 to mass determination, see [12-14] 
* Typically the signal has a slightly larg er p^ w 
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Figure 1: The distributions of M?2 for various values of run- 



As the transverse mass is bounded above by the invariant mass, the above relation 
implies 



M T2 < min(M (1) ,M (2) ) < 



(7) 



where M' 1 ''' 2 ' are the invariant masses of the intermediate W bosons. One then immedi- 
ately finds that the Mt2 of dileptonic Higgs decay is bounded as (upon ignoring the finite 
decay widths of the Higgs and W bosons) 



M T2 < mm(M w 



ma- 



in Fig. [Tl we depict the Mt2 distribution for several different values oimn, where the event 
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set is generated for the process gg — > H — > W^W 7 — > ZvZV at the LHC beam conditiorj§. 
The results show that the endpoint value is indeed m#/2 for m# < 2Mw, while it is Mw 
for m# > 2Myy. 

The MAOS momenta [7] are the four-vector variables which approximate the momenta 
of invisible particles in collider event of the type 

X(p + k)X'(q + I) - 7(p) + *(*;) + V'(q) + X '(l), 

where X and X' denote pair-produced (not necessarily identical and not necessarily in on- 
shell) mother particles, V(p) and V'(q) represent (one or more) visible particles with total 
momenta p and q, which are produced by the decays of X and X', respectively, and x an d 
x' are two invisible particles having the same mass. The transverse MAOS momenta are 
defined as the trial transverse momenta of x an d x' that determine the Mti of the event, 
i.e. k/p and It minimizing max {Mt(X), Mt(X')} under the constraint kr + lr = \/>t, where 
Mt(X) and Mt(X') are the transverse masses of X and X', respectively. For the process 
WW —> l(p)v{k)V (q)v' {I) , the transverse MAOS momenta are uniquely determined by the 
conditions 



M T2 (p T ,q T ,pV) = ^2(|p T ||k^ s | -p T -k^ aos ) 



9 (\ n 111 maos I „ 1 maos \ 
„J i,maos i lmaos (ci\ 

p T = k r + 1 T . (9) 

For an event without initial state radiation, i.e. an event with \i T 7W = 0, we have pV — 
— (Pt + Qt) and M?2 given by (jSJ). In this case, the transverse MAOS momenta are simply 
given by 

k™ aos = -q T , l™ aos = -p T . (10) 

There can be two different schemes to define the longitudinal MAOS momenta. One 
is to require the on-shell conditions for both the invisible particles in the final state and 
the mother particles in the intermediate state. In the case of WW — > l(p)v(k)V (q)v' (I) , it 
corresponds to 

(^rnaos) = (^maos) = 0, (j) + A; maos ) = (q + / ma os) = (I I) 

which results in 



IPtT 

/maos / I \ 

L { > ~ T7T12 
Or 



p L A±J\p T \2+pl J A 2 - \p T \ 2 \k^ 



,maos 1 2 



q L B± J\q T \ 2 + q\ J B 2 - |q T | 2 |l? aos | 2 



(12) 



§For the analysis in this section, we do not include the effects of hadronization and detector smearing, 
while those effects are incorporated in the analysis of the next section. 
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where A = M w /2 + p T • k^ aos and B = M w /2 + q T ■ l?J aos . Note that the above longitudinal 
MAOS momenta have four-fold degeneracy for each event (two-fold degeneracy for each 
neutrino MAOS momentum). Another possible scheme is to require 

(fcmaos) = (^maos) = 0, (j? + A; maos ) = (q + / ma os) = Mrp 2 , (13) 

which gives unique longitudinal MAOS momenta as 

1 1 maos I 

ir = -t^ql. (14) 

To distinguish these two schemes, the MAOS momenta of (fTUj) and (fT2l will be called the 
original MAOS momenta as they are the one originally defined in [7], while ({TO]) and (fl4l) 
will be called the modified MAOS momenta. 

A nice feature of the MAOS momenta is that they provide a systematic approximation 
to the invisible particle momenta [7], i.e. the neutrino momenta in our case. For the case 
of rriH > 2Mw, both of the two W bosons in H — *• W + W~ are on mass-shell, and then the 
endpoint value of Mt2 is given by Mw- In this case, one easily finds that both the original 
MAOS momenta and the modified MAOS momenta approach to the true neutrino momenta 
in the limit of the endpoint event with M T2 = M w . (Note that kf aos (+) = kf aos (-) for the 
endpoint event.) For generic events with Mt2 < Mw, the MAOS momenta generically differ 
from the true neutrino momenta. Even in these cases, we can infer from its distribution that 
the MAOS momentum provide a reasonable approximation to the true neutrino momentum. 
As the approximation gets better for larger value of Mt2, one can systematically improve 
the efficiency of approximation with Mt2 cut. 

In Fig. [2J we show the distributions of the difference between the MAOS momentum 
and the true neutrino momentum, 

At — h trne iLmaos 
AKt,L — K T,L ~ T,L ) 

for tuh = 180 GeV and m# = 140 GeV, respectively. The left panels include the distribu- 
tions of the full event set for gg — > H — > WW — > Ivl'v 1 generated at the LHC condition, 
while the right panels show the distributions of the top 10% subset near the endpoint 
of M?2- By definition, the original and modified MAOS schemes give the same trans- 
verse MAOS momenta, thus the same Akr distribution (black). For A/c^ in the original 
MAOS scheme, we construct its distribution using the two solutions {fc™ aos (+), Z™ aos (+)} 
and {/c™ aos (— ), Z™ aos (— )} for each event. For m# > 2Myy, our results indicate that the 
original scheme (red) is a bit better than the modified scheme (blue), if one considers the 
full event set. However, the two schemes show similar performance if one employs a proper 
M T2 cut, say a cut selecting about 30% of the near-endpoint events. Both the modified 
and original schemes recover the true neutrino momentum for the exact endpoint event 
with Mt2 = Mw- On the other hand, when ran < 2Mw so that one or both VT-bosons 



l,maos 
r,maos _ IT I 

Pt 
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Figure 2: Distributions of Ak? = k^ ue — fc™ aos (black), Ak^ for the original MAOS momenta 
(red), and Aki for the modified MAOS momenta (blue). The upper frames are for m# = 
180 GeV, and the lower frames are for tuh = 140 GeV. Full event set is used for the left 
frames, while the top 10 % near endpoint events of Mt 2 are used for the right frames. 

are in off-shell, the modified MAOS is clearly the better choice to approximate the neu- 
trino momenta. For off-shell W, the original MAOS scheme does not give correct neutrino 
momenta even for the endpoint event of Mt2, while the modified MAOS scheme does. 

Once the MAOS momenta of neutrinos are obtained, one can construct the MAOS 
Higgs mass: 

(m- aos ) 2 = (p + fc maos + q + / maos ) 2 . (15) 

A nice feature of m^ aos is that its distribution has a peak at the true Higgs boson mass, 
which becomes narrower under a stronger Mt2 cut. In Fig. [3l we show the MAOS Higgs 
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Figure 3: The MAOS Higgs mass distributions for the full event set (left frames) and the 
top 10 % near endpoint events of Mt2 (right frames) in the original (red) and modified 
(blue) MAOS schemes. The upper frames are for run = 180 GeV and the lower frames are 
form H = 140 GeV. 

mass distributions of the full event set (left panel) oigg — > H — > WW — > Ivl'v' generated for 
ma = 180 and 140 GeV at the LHC condition, and also of the top 10% near endpoint events 
of Mt2 (right panel), for both the original MAOS scheme (red) and the modified MAOS 
scheme (blue). Our results suggest that the modified MAOS scheme can provide a good 
approximation to the invisible neutrino momenta for both m# < 2Mw and ra# > 2Mw, 
under a suitable M?2 cut selecting near endpoint events. In fact, such an M^ 2 cut is useful in 
another sense as it enhances the signal to background ratio. As is well known, the dileptons 
from the Higgs decay are likely to have smaller opening angle than the background, which 
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is essentially due to that the Higgs boson is a spin zero particle. Then, the expression ([6]) 
of Mt2 indicates that the dileptons from the Higgs decay are likely to have larger Mt2 than 
the background because of the same reason. 

In the next section, we will perform the likelihood analysis for the MAOS Higgs mass 
distribution to determine the true Higgs boson mass, while including some of the main 
backgrounds as well as the detector effects. We will use the modified MAOS scheme since 
the original MAOS scheme does not work for the case of m# < 2Mw, while the modified 
MAOS scheme works well for both mn < 2Mw and m# > 2Mw- 

3 MAOS reconstruction of the Higgs boson mass 

To investigate the experimental performance of the MAOS reconstruction of the Higgs 
boson mass at the LHC, we use the PYTHIA6 .4 Monte Carlo (MC) generator at NLO [15]. 
The generated events have been further processed through the fast detector simulation 
program PGS4 [16] to incorporate the detector effects with reasonable efficiencies and fake 
rates [17]. Assuming the integrated luminosity of 10 fb _1 , we have generated the MC event 
samples of the SM Higgs boson signal and the two main backgrounds. For the signal, we 
consider the Higgs boson production via the gluon fusion: gg — > H . For the Higgs mass 
range 130 GeV < mn ^ 200 GeV, the produced Higgs boson decays mainly into a pair 
of W bosons. We take into account all the dileptonic decay channels of the W bosons 
to enhance the signal, W — > Iv with I — e, \i. The dominant background comes from the 
continuum qq , gg — > WW — * Ivl'v' process, and we include also the tt background in which 
the two top quarks decay into a pair of W bosons and two b jets. 

Following [1], we have imposed the following basic selection cuts on the Higgs signal 
and the backgrounds: 

• Require that the event has exactly two isolated, opposite-sign leptons (electron or 
muon) with > 15 GeV and \rj\ < 2.5. 

• 12 GeV < mu < 300 GeV. 

• \$ T \ > 30 GeV. 

• No b jets. 

• No jets with p T > 20 GeV. 

It is well known that the background can be significantly reduced by exploiting the 
helicity correlation between the charged lepton and its mother W boson. Introducing the 
transverse opening angle between two charged leptons, A <!>/;, the Higgs signal tends to have 
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a smaller A$;/ than the background, which is essentially due to the fact that the Higgs 
boson is a spin zero particle. Selecting the events with large value of Mt2 similarly enhances 
the signal to background ratio, which can be understood by the correlation between My 2 
and in ([6]). In our case, this M^ 2 cut is particularly useful since it enhances also the 
accuracy of the MAOS reconstruction of the neutrino momenta as discussed in the previous 
section. In Fig. HI we show the scatter plots of Mt2 and A§u for the signal (left panel) 
and the background (right panel), obtained after imposing the above basic selection cuts 
to the data set for m H = 170 GeV, while including the detector effects. As anticipated, 
the low A$;; and high M^ 2 region is more populated by the signal events, while the high 
A&U and low Mt2 region by the backgrounds. We can also notice a correlation between 
A&u and Mt2 suggested by (EJ). Fig. [5] shows the Mr 2 distribution (left panel), again for 
rrin = 170 GeV, of the events with A$;/ < 1.6, and the A$« distribution (right panel) 
of the events with M^ 2 > 67 GeV, where the shaded regions represent the backgrounds. 
Here, the tail of the Mp 2 distribution beyond My/ is mainly due to the VF-boson width. 
We observe that the signal is more likely to have larger M^ 2 . Furthermore, the right panel 
of Fig. [5] indicates that the Mt2 cut can significantly enhance the efficiency of the A$^ cut. 
We thus introduce two additional cuts: 

• A$« < A$f\ 

• M T2 > M$, 

where A$^ ut and M^ 2 * are chosen to optimize the Higgs mass measurement using the MAOS 
mass distribution, and their values for various m H are listed in Table [TJ 



Table 1: The A$^ and My 2 cuts for various values ofmu- 



m H (GeV) 


130 


140 


150 


160 


170 


180 


190 


200 


A$r 

Mj$ (GeV) 


1.85 
38.0 


1.70 
51.0 


1.65 
57.0 


1.50 
66.0 


1.60 
67.0 


1.70 
68.0 


1.90 
69.5 


2.05 
70.0 



In Table [21 taking the case of m# = 170 GeV as a specific example, we show how 
the numbers of signal events and background events are changing under each selection 
cut. Comparing with the ATLAS cut flows reported in [1], the signal and the dominant 
WW background are in excellent agreement except for that the overall number of events is 
somewhat larger in our cas^l This may be attributed to the differences in simulating the 

^Note that we consider all of the ee, /x/x, e/z events, while [1] included only the e\i events. 
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Figure 4: The scatter plots of Mj- 2 and A<&u for the signal (left) and background (right) 
events after imposing the basic selection cuts. For the signal, tuh = 170 GeV is taken. 




Figure 5: The My 2 distribution (left) for the events with < 1.6, and the A$» distri- 
bution (right) for the events with Mt2 > 67 GeV. Signal events are taken for rrin = 170 
GeV, and the shaded regions represent the backgrounds. 

jet reconstruction, detector effects and triggering. For the ti background, we have a sizable 
number of events even after imposing the basic selection cuts, but it can be suppressed to 
a negligible level by the A$« and Mt2 cuts. 

In Fig. El we show the m^ aos distributions for various (nominal) values of the Higgs 
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Table 2: Cut flows for m H = 170 GeV with A$Jf = 1.6 and M^f = 67 GeV at 10 fb' 1 . 



Selection 


Selection cuts 


99 


WW 


tt 




TiPrtfon ^plppfion -1-777 77 


A AA5 


18 501 


139 256 


Basic 


\$ T \ > 30 GeV 


4,012 


12,801 


120,597 


Selection 


b- veto 


3,956 


12,656 


60,438 




Jet veto 


2,039 


8,096 


1,287 


Tuned 
Selection 


A$„ < A<^ ut 

M T2 > 


1,621 
619 


2,939 
585 


332 
107 



boson mass, obtained by the modified MAOS momenta of neutrinos, while incorporating 
the Mt2 and A&u cuts listed in Table [H We observe that the modified MAOS scheme is 
working nicely and each distribution has a clear peak at the true (nominal) Higgs mass 
independently of whether m# is below or above 2Mw- 

With the m^ aos distribution constructed as above, we performed a template fitting 
to determine the Higgs boson mass. Here a template means a simulated distribution with 
a trial Higgs mass which, in general, is different from the nominal one used to generate 
the data. For each distribution with the nominal Higgs mass ma, the 11 templates are 
generated with the trial Higgs masses between tuh — 10 GeV and rrin + 10 GeV, in steps 
of 2 GeV. For example, in each frame of Fig. [61 the solid line shows the template when 
the trial Higgs mass is the same as the nominal one. Each template is normalized to the 
corresponding nominal distribution. 

The likelihood between a nominal data distribution and a template is defined as the 
product of individual Poisson probabilities computed in each bin i over the M bins in the 
fit range: 

„ -^Ce-^m™ 1 . , 

^n^T^. (is) 

where rij and rrii denote the number of events in the z-th bin of the nominal distribution and 
the normalized template, respectively. In Fig. [TJ we show the log likelihood distribution for 
various Higgs masses in the range between 130 GeV and 200 GeV. The solid line shows the 
result of a quadratic fitting for each value of mjj. The fitted Higgs boson masses together 
with 1-cr error are listed in Table [3] (see also Fig. [8]) for various input Higgs masses, where 
the l-o" deviated value is defined as the one increasing — ln£ by 1/2 [18]. 
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Table 3: The fitted Higgs boson mass for various input values of ran at 10 fb 1 . 



M H (GeV) 


130 


140 


150 


160 


170 


180 


190 


200 


Fitted value (GeV) 
l-o" error (GeV) 


130.0 
2.4 


140.1 
1.7 


150.9 
1.2 


160.6 
1.0 


170.3 
0.9 


179.4 
1.4 


190.4 
2.0 


199.7 
3.5 



4 Conclusions 

In this paper, we have examined the prospect of measuring the Higgs boson mass using 
the MAOS neutrino momenta in the Higgs decay H — > WW — > lulu. To optimize the 
efficiency of the method, we have employed an event selection to combine the well known 
dilepton azimuthal angle cut, < A$™*, with an Mj- 2 cut selecting only the events with 
Mt2 > Mt£. This Mt2 cut enhances the efficiency of the MAOS momentum approximation 
to the true neutrino momentum, and the signal to background ratio also. Under such 
selection, the MAOS Higgs mass distribution constructed with the measured charged lepton 
momenta and the MAOS neutrino momenta shows a clear peak at the true Higgs boson 
mass. Likelihood fit analysis for the MAOS Higgs mass distribution suggests that it can 
provide a precise determination of the Higgs boson mass, and also be useful for the discovery 
or exclusion of the Higgs boson in certain mass range. 

Our analysis can be improved in many respects, e.g. with more extensive study of 
backgrounds, with more complete detector simulation, and with multi-dimensional fitting 
including other observables. At the moment, it is not straightforward to compare our 
results with those of Ref. [1] providing a more complete analysis using the transverse mass 
variable M^ pprox = M T {m vv = m u ) (together with A$ u and pf w ) instead of the MAOS 
Higgs mass m# aos , where Mf = m\ + m 2 vl/ + 2(E^E^ U — p" • p™) denotes the transverse 
mass of WW — > lulu. On the other hand, it is relatively easy to compare our analysis 
with [10] which discusses the Higgs mass determination with M^ me = Mt{v(i vi , = 0) in a 
simple context including only the dominant background qq — > WW without taking into 
account the detector effects. Our results appears to be better than (or comparable to) 
those of [10]. This indicates that incorporating the collider variable m^ aos might improve 
significantly the accuracy of the Higgs mass measurement, as well as the significance of the 
Higgs boson discovery or exclusion. 
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Figure 6: The m^ aos distribution for the nominal data (dots) and the template (solid line) 
for various Higgs boson masses: van = 140, 150, 160, 170, 180, 190 GeV from the top left 
to the bottom right. In each frame, shaded region represents the backgrounds. 
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Figure 7: The relative log likelihood distributions for various Higgs boson masses. The solid 
line shows the result of the quadratic fitting for each value ofm H . e maximum 

likelihood which was determined as the minimum of a fit to the — In C distribution. 
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Figure 8: The band showing the 1-a deviation error for the Higgs boson mass determined 
by the m^ aos distribution. The dots and lines denote the Higgs boson mass obtained by the 
likelihood fit. 
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